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Protecting and Enhancing Spin Squeezing via Continuous Dynamical Decoupling 
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Realizing useful quantum operations with high fidelity is a two-task quantum control problem 
wherein decoherence is to be suppressed and desired unitary evolution is to be executed. The 
dynamical decoupling (DD) approach to decoherence suppression has been fruitful but synthesizing 
DD fields with certain quantum control fields may be experimentally demanding. In the context 
of spin squeezing, here we explore an unforeseen possibility that continuous DD fields may serve 
dual purposes at once. In particular, it is shown that a rather simple configuration of DD fields 
can suppress collective decoherence and yield a 1/N scaling of the squeezing performance (A'^ is the 
number of spins), thus making spin squeezing more robust to noise and much closer to the so-called 
Heisenberg limit. The theoretical predictions should be within the reach of current spin squeezing 
experiments. 
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The feasibility of "dynamical decoupling" (DD) [l[ in 
effectively isolating a quantum system from its environ- 
ment has attracted great theoretical and experimental in- 
terests. Remarkable progress towards efRcient protection 
of quantum states has been achieved [2|-l5[. In contrast, 
high-fidelity protection of quantum operations (such as 
quantum gates or quantum metrology schemes) is more 
challenging experimentally. As decoherence must be sup- 
pressed during quantum operations, it is natural to syn- 
thesize DD fields with other fields implementing a de- 
sired quantum operation. However, this bottom-up ap- 
proach may require complicated coherent control fields. 
For example, explicit solutions to dynamically corrected 
quantum gates arc sophisticated [6|-|8| , and even a simple 
quantum metrology protocol, when combined with DD, 
already becomes a rather involving practice |9| . 

A top-down approach to the protection of useful quan- 
tum operations should be a worthy direction, along which 
we aim to better exploit the system's own Hamiltonian 
under DD fields. In essence we are faced with a two-task 
problem: decoherence is to be suppressed and desired 
(almost) unitary evolution is to be executed. Is it possi- 
ble to directly construct DD fields serving the dual tasks 
at once? Motivated by recent studies of decoherence ef- 



feet 



s on spin squeezing 



^ and by recent exciting 

experiments of spin squeezing |18l . |l9[ , here we use the 
spin squeezing context to give a positive answer to our 
question. That is, it is feasible to protect and enhance 
spin squeezing at the same time by searching for a special 
configuration of DD fields. The dual roles of DD fields 
arise from two facts: (i) DD fields modulate both system- 
bath interaction and the system Hamiltonian itself that 
describes the spin-spin interaction, and (ii) the system 
Hamiltonian itself under the modulation of DD fields may 
generate more useful quantum evolution. We empha- 
size that the enhancement in spin squeezing we achieve 
is not a secondary outcome of decoherence suppression. 



Rather, the enhanced spin squeezing is far superior to 
what can be normally achieved under "decoherence-free" 
conditions. Indeed, we predict a 1/N scaling of the ob- 
tained spin squeezing performance, where N is the num- 
ber of spins. Our theoretical results should be testable 
by modifying existing experiments. In addition, because 
spin squeezing is closely related to multi-partite entan- 
glement [201, l21| , the results are also of interest to ongoing 
studies of entanglement protection 22l426l| . 



Consider then a collection of N identical spins (or 
qubits). Ill terms of the standard Pauli matrices, the 
dynamics can be described by collective angular momen- 
tum operators Jk = ^^m=i"k"^ with k = x,y,z. In 
spin squeezed states, quantum fluctuations of the collec- 
tive angular momentum in one direction arc significantly 
reduced at the price of increased uncertainty in another 
direction 271 1 2 Si. th us offering higher precision in quan- 
tum metrology j29l |30[. For instance, we can improve 
high-precision spectroscopy and atomic clocks which are 



currently limited by spin noise [31 



Using the angular momentum commutation relations, 
the two most widely used measures of spin squeezing ^^ 
23 and ^1 [ill, are found to be bounded by 1/7V 



This fundamental limit to the amount of spin squeezing 
achievable reflects the Heisenberg precision limit in quan- 
tum measurement |28l. l33l|. In practice, the achievable 
degree of squeezing is considerably worse than the 1/iV 
limit for two main reasons. First, squeezing is in general 
degraded by decoherence or noise [3J] . The environment 
tends to destroy squeezing, causing the sudden death of 
squeezing [l3|. It may also change the optimal squeez- 
ing time window in an unpredictable way, leading to a 
non-optimal squeezing generation [19| . Second, Hamilto- 
nians that can be implemented so far cannot reach the 
1/7V scaling in theory. For instance, in two recent experi- 
ments [l8|, ll9J based on two- mode Bose-Einstein conden- 
sate (BEC), the so-called one-axis twisting (OAT) Hamil- 



tonian [27[ Hqat = X^x i^ realized, which can at most 
generate £,j^ ~ I/N3 , not to mention decoherence effects. 
It is thus clear that protecting spin squeezing against de- 
coherence and pushing spin squeezing towards the 1/A^ 
scaling would be of wide interest. 

We start by considering a OAT system interacting with 
an environment, with the total Hamiltonian modeled by 
H = Hq + Hb + ffsB, where, Hq = Hqat = xJl, Hb is 
the Hamiltonian of the environment, and -ffsB represents 
the system-environment coupling, -ffse is assumed to be 



H^ 



SB 



^x^ X 



Byjy 



BzJz, 



(1) 



where the Bk are arbitrary bath operators (or randomly 
fluctuating noise for a classical bath). Though coupling 
terms that are nonlinear in Jk are not considered here, 
the HsB in Eq- (HI is already quite general insofar as it 
covers a broad class of problems with both dephasing and 
relaxation. 

There is a standard route to seek a control Hamilto- 
nian Hc{t) that can effectively average out Hsb and hence 
suppress decoherence. In particular we consider a contin- 
uous Hc{t) of period tc-, whose time-ordered exponential 
defines a unitary operator Uc{t) — Texp[— i/^ Hc{t')dt'] 
(fi = 1 throughout), with Uc{t+tc) ^ Uc{t). The Magnus 
expansion [35|, |36[ indicates that if 



Ul{t)HsBUcXt)dt^Q, 



(2) 



then to its first order H^b is suppressed. Extending 
previous studies for singlc-qubit and two-qubit systems 



371 138| , we choose 

Tj t±\ __ -2-ninyJyt/tc —2iTin^J^t/tc loX 



where Ux and Uy are non-zero integers. For any Ux ^ ny^ 
Uc{t) in Eq. ([3|) satisfies the first-order DD condition of 
Eq. ([2]). Qualitatively, such Uc{t) causes the collective 
angular momentum operators to rapidly rotate in two 
independent directions and as a result, ifsB is averaged 
out to zero. Using idUc{t)/dt = Hc{t)Uc{t), we obtain 
the following DD control Hamiltonian, 

Hc{t) ~ LUHyJy + Lunx[Jx cos{Lunyt) — Jz sm{Lonyt)], (4) 

where ui = lix jtc- The total system Hamiltonian Hs{t) ^ 
Hqat + He (i) then becomes 

Hs{t) = xJx +^nyJy+^nx[Jx sin(ajrij,t) — J^ cos{ujnyt)]. 

To elaborate how Hgit) can be realized, we rotate the 
coordinate system along the y-axis by 7r/2, transform- 
ing Hs{t) to H'g{t) = xJz +ojnyJy — ujnx[JxSm{,<^nyt) + 
Jz cos{ujnyt)]. We now comment on each term of H'g{t). 
The first Jj term describes spin-spin interaction, as is 
realized in experiments |18l . I19j . The last term linear in 
Jz can be realized by an oscillating energy bias using for 



example a time-dependent Zeeman shift. The Jx and Jy 
terms can be generated by use of electric-dipole inter- 
action - considering a circularly polarized transition, a 
constant electric field along y direction and an oscillat- 
ing field along x direction lead to the desired Jx and Jy 
terms [H. 

With a continuous control Hamiltonian Hc{t) imple- 
mented, decoherence can be well suppressed for suffi- 
ciently large w. Two observations are in order. First, 
as shown in Eq. (|4]), infinite DD solutions with differ- 
ent (nxjUy) combinations are found. Second, the con- 
trol Hamiltonian averages out Hq,b via fast modulations 
of Jfc, so the system's self-interaction term Jj is nec- 
essarily modulated at the same time. One opportunity 
is then emerging: among all the DD solutions, can we 
identify a particular type that modulate the system's 
self-interaction Hamiltonian in a useful manner so as to 
enhance squeezing while suppressing decoherence? 

With the system decoupled from its environment, it 
can be shown that the system evolution operator is given 
by Us{t) ~ Uc(t)e~'^^* [36|, where the time-averaged 
Hamiltonian H is found to be 



^-f I ^ Ul{t)JlU,{t) dt. 

'c Jo 



(5) 



A straightforward though rather tedious calculation 
yields 

Ul{t)JlUS) 
1 



&va.{2ijjnxt) sin {ujnyt)[JzJy + JyJz 



H sm.{2ujnyt) coii{u!nxt)[JxJz + JzJx] 

-\ sin(2a;n.yt) sa\{ujnxt)[JxJy + JyJx] 

+ jj cos^ {ujUyt) + Jy sin [ujUxt) sin {uiriyt) 



J^ cos (ujUxt) sin {ujnyt). 



(6) 



Using Eq. (jH]), one finds the time-averaged Hamiltonian 
H has two different forms. Specifically, if Ux 7^ Srij,, 
H = jJ^ (up to a constant), which is just the original 
OAT Hamiltonian with the nonlinear coefficient scaled 
down by a factor of four. If Ux = 2ny, which we call the 
"double-resonance" (DR) condition, we obtain (up to a 
constant) 



H = Hi 



DR. 



fw 



'Jx'Jy ~r JyJxj 



(7) 



Remarkably, i?DR is seen to be a mixture of a OAT 
Hamiltonian and a well-known two-axis twisting (TAT) 
Hamiltonian iJxAT = xiJxJy + JyJx) [271, |28|. Since 
Htat is known to produce the best scaling of squeezing, 
we are motivated to examine the squeezing performance 
of Hbr, naturally obtained by one type of DD fields to 
fight against both relaxation and dephasing. 
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FIG. 1. (color online) Spin squeezing measure ^| against time 
t using the OAT Hamiltonian (dot-dashed, red), Hun (solid, 
dark blue), and Hs{t) with A'^cyc ~ 5 (dotted, magenta) and 
Ncyc = 20 (dashed, blue) for iV = 10 (i.e., J = 5). Note 
that the dynamics generated by Hs{t) with A'^cyc = 20 are al- 
most indistinguishable from the dynamics generated by Hub.- 
Here we use Ux — 2 and Uy = 1, and imin was found to be 
approximately 0.491. 



FIG. 3. (color online) For J — 5, spin squeezing measure ^| 
against time t using the bare OAT Hamiltonian without noise 
(dot-dashed, red), Hok (solid, dark blue), the OAT Hamil- 
tonian with noise but without DD fields (dashed, magenta), 
and the OAT Hamiltonian in the presence noise and the DD 
fields with n^ — 2, Uy — 1, imin ~ 0.491, and A'cyc = 20 (dot- 
ted blue line, which is almost on top of the solid line). An 
average over 2000 sample paths of the noise was taken. The 
noise parameters are a = 2 and a^ = 20. 



We consider an initial state \J,—J) describing all spins 
"pointing dovi^n" , that is, an eigenstate of J^ with eigen- 
value ^N/2. To verify our expressions of /?dr and in- 
vestigate its potential benefits we first switch off H^b 
and evolve the wavefunction numerically under Hs{t) = 



HoAT + Hc{t). We used 



_ 4min(AJH^)^ 

N 



[27| to quantify 



squeezing, where fi± denotes a direction perpendicular to 
the mean spin direction and the minimum is taken over 
all such directions. The dynamical behavior of ^ is found 
to be essentially the same as ^ (^ ~ ^s{J/\{J)\)'^)! so 
only the behavior of ^ is presented below. We set x = 1 



and tc ~ in 



i/iVcyc, where imin 



is the optimal squeezing 




0.12 



FIG. 2. (color online) Same as in Fig. [T] but this time we 
have N = 100 (i.e., J = 50) with iVcyc = 10 (dotted, ma- 
genta), Ai'cyc = 30 (dashed, blue), and imin ~ 0.0909. It 
is also observed here that the dynamics generated by Hs{t) 
with A'cyc = 30 are well captured by the dynamics generated 
by Hob. (solid, dark blue). 



time for Bur, and A'cyc is a positive integer. In Fig.[TJ the 
squeezing performance of H^ (i) is compared with that of 
Hqat = xJx^ for iV = 10. It is seen that Hs{t) generates 
much better squeezing - the minimum value of ^ using 
Hs(t) is approximately 0.15, but we can only achieve a 
value of approximately 0.2 using Hqat- In addition, it is 
seen from Fig. [T] that the dynamics under Hs{t) is indeed 
well captured by the dynamics under the time-averaged 
Hamiltonian i?DR- In particular, within each period of 
tc, Hs{t) and H^n yield some fluctuating differences in 
£,g, but at integer multiples of tc when Uc{t) = 1, excel- 
lent agreement between them is observed. Further, if we 
reduce tc by increasing the strength and the frequency of 
the periodic control Hamiltonian Hc{t), the fluctuating 
differences become smaller. This reflects the fact that 
iJoR is obtained under a first-order approximation. We 
next increase A^ tenfold and essentially the same results 
are obtained in Fig. ^ but with the advantage gained by 
our DD fields displayed even more evidently. For exam- 
ple, in Fig. 2 it is seen that the minimum values of ^ us- 
ing Hs{t) and Hqat arc approximately 0.019 and 0.048, 
respectively - so more than a two-fold improvement can 
therefore be obtained. Other calculations indicate that 
for larger A^, the accuracy /?dr goes down with fixed 
tc- Thus, a higher driving frequency a; would be more 
favored as A^ increases. 

Having confirmed that the double-resonance condition 
rix = 2ny is useful for spin squeezing, let us now turn to 
the full problem by switching on the system-environment 
coupling. For convenience we model B^, By and B^ in 
Eq. ([T]) as three independent Gaussian colored noise pro- 
cesses, with the same inverse correlation time a and noise 



variance a^ . Wc numerically compute the dynamics of 
squeezing for Hs(t) in the presence of noise and then 
compare it with that generated by ^oatj with noise or 
without noise. As shown in Figs. [3] and 21 Hs{t) with 
noise yields much better squeezing than Hqat with noise. 
This may be understood as an outcome of decoherence 
suppression. On the other hand, Hs{t) with noise also 
generates better squeezing than Hqat in the absence of 
noise. Hence our DD fields have played one more role in 
addition to decoherence suppression. Note also that the 
time to obtain maximum squeezing is in excellent agree- 
ment with that obtained from -ffDR, hence avoiding de- 
coherence effects on the optimal squeezing time and also 
confirming again the usefulness of i^DR in predicting the 
optimal squeezing time. 

Having shown how DD fields may suppress decoher- 
ence and enhance the spin squeezing generation as a uni- 
tary process, we finally investigate how the squeezing 
performance of Hbr scales with N. Within the valid- 
ity regime of Hur as an effective Hamiltonian (for de- 
scribing the dynamics associated with the OAT Hamilto- 
nian in the presence of noise and continuous DD fields) , 
the scaling of the squeezing performance of -ffDR with 
N represents to what degree our DD fields can protect 
and enhance spin squeezing. Calculations for even larger 
values of N then become necessary. We first compare 
the performance of 7?dr with what is known to give 
the best scaling behavior, namely, the TAT Hamiltonian 
Htat = xiJxJy + JyJx)- Significantly, although i^DR 
produces slightly less squeezing than Htat, two close 
and parallel lines describing their respective performance 
are seen in Fig. [S] indicating that both cases give the 
^1 ~ 1/A^ scaling. By contrast. Fig. [5] also presents the 
l/N^/^ scaling of the OAT Hamiltonian [also see the inset 
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FIG. 4. (color online) Same as in Fig. (3] but now with J = 50, 
imin ~ 0.0909, A'cyc = 30, and a^ = 100. An average over 
100 sample paths of the noise was taken. Similar to what 
is observed in Fig. [31 the squeezing performance of the OAT 
Hamiltonian in the presence of noise and DD fields (dotted, 
blue) is much better than that of a bare OAT Hamiltonian in 
the absence of noise (dot-dashed, red) . 



of Fig. [5] for a comparison of two different scalings] . The 
DD fields under the double-resonance condition hence al- 
lows squeezing to occur in the presence of noise and in the 
mean time brings about a squeezing enhancement factor 
of N^l'^ , which is in principle unlimited as N increases. 

It is also interesting to note how this work differs from a 
recent proposal for realizing TAT Hamiltonian by apply- 
ing a designed pulse sequence to a OAT Hamiltonian |40[ . 
While our starting point is continuous DD fields for deco- 
herence suppression, the short control pulses considered 
in Ref. |40| do not average out i^ss in Eq. (1) to zero. 
Further, the effective Hamiltonian iJoR found here un- 
der a double-resonance condition is a mixture of OAT 
and TAT Hamiltonians. To our knowledge, i^DR is a 
newly found, physically motivated Hamiltonian that can 
1/A^ scaling. 
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FIG. 5. (color onhne) Minimum value of ^g - in dB, defined as 
lOlog(^l) - plotted against N (log scale) for dynamics gener- 
ated by the OAT Hamiltonian (upper, dot-dashed, red), the 
TAT Hamiltonian (lower, dashed, magenta), and the time- 
averaged Hamiltonian Hdr (middle, solid, green). In the in- 
set, the upper (dot-dashed, red) line has been transported 
vertically to show clearly the better scaling behavior of -ffDR 
and the TAT Hamiltonian. The behavior of ^^ (not shown) 
is similar. 



To conclude, by considering a class of continuous fields 
to suppress both dephasing and relaxation in the dynam- 
ics of spin squeezing, we are able to identify a special type 
of DD solutions that can effectively yield a previously un- 
known spin squeezing Hamiltonian, generating the 1/iV 
scaling of squeezing performance in the presence of an 
environment. With their dual roles in decoherence sup- 
pression and in generating more useful quantum evolu- 
tion identified, the found DD fields are appealing from an 
experimental point of view. Our results should be able 
to help design new experimental studies of spin squeezing 
based on one-axis twisting Hamiltonians (such as those 
using two- mode BEG). Indeed, by exploiting system's 
own spin-spin interaction Hamiltonian under the mod- 
ulation of continuous DD fields, we expect to see other 



interesting DD designs that can carry out desired quan- 
tum operations while protecting quantum coherence. 

A.Z.C. would like to thank Derek Ho and Juzar 
Thingna for helpful discussions. 
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